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DIRECTED POLYMERS AND THE QUANTUM TODA 

LATTICE 

By Neil O'Connell 
University of Warwick 

We characterise the law of the partition function of a Brownian 
directed polymer model in terms of a diffusion process associated with 
the quantum Toda lattice. The proof is via a multi-dimensional gener- 
alisation of a theorem of Matsumoto and Yor concerning exponential 
functionals of Brownian motion. It is based on a mapping which can 
be regarded as a geometric variant of the RSK correspondence. 

1. Introduction. Let Bi{t), B2{t), . . . , B]\f{t),t > be a collection of 
independent standard one-dimensional Brownian motions and write Bi{s, t) = 
Bi{t) — Bi{s) for s <t. Let /3 G R, t > 0, and consider the random variable 

(/5) = / e/3(Bi(.0+B2(.i,.2)+-+Biv(.iV-i,t))^g^ . . . 

Jo<si<-<SN-i<t 

This is the partition function for a model of a 1 -|- 1 dimensional directed 
polymer in a random environment which has been introduced and studied 
in the papers [41, 42, 53]. The free energy density is given explicitly by 

lim 1 logZ^(/3) = mf Wh - ^{t)] - log/S^, 

N^oo iV t>0 

almost surely, where "^{z) = r'{z)/T(z). The law of Z/^(/3) is well-understood 
in the zero temperature limit /3 — t- oo, where it has close connections with 
random matrices. Define 

= lim ilogZf(/3) 

(1) = max {Bi{si)+B2{si,S2) + --- + BN{sN~i,t)) 

0<si<---<sjv-i<t 

Note that, by Brownian scaling, the law of /\/t is independent of t. 

Theorem 1.1 The random variable has the same distribution as the 
largest eigenvalue of a NxN GUE random matrix [2, 23]. In fact [7, 43], the 
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stochastic process {Mf, t > 0) has the same law as the largest eigenvalue 
of a standard Hermitian Brownian motion, that is, it has the same law 
as the first coordinate of a Brownian motion conditioned (in the sense of 
Doob) never to exit the Weyl chamber Cn = {x E : xi > • • • > x^}, 
started from the origin. This is a diffusion process in C n with infinitesimal 
generator A/2 + V log h ■ V where 

(2) h{x)= Yl 

This connection with random matrices yields very precise information con- 
cerning the distribution and asymptotic behavior of when N is large. 
For example, it follows that 

lim P (m^ <2N + xN^/'^\ = F2{x), 

where F2 is the Tracy- Widom distribution [56]. 

In this paper we obtain an analogue of Theorem 1.1 for the stochastic 
process (logZj^(/3), i > 0). We will show that, for each /5 > 0, this process 
has the same law as the first coordinate of a diffusion process in which 
is closely related to the quantum Toda lattice. This yields an analytic de- 
scription of the law of Zl^{(3) which should provide a good starting point 
for further asymptotic analysis. 

2. The quantum Toda lattice. The quantum Toda lattice is a quan- 
tum integrable system with Hamitonian given by the Schrodinger operator 

(3) ^ = Ea^-2E^""^""- 

i=l * i=l 

It is closely related to the Lie algebra g/^: the exponents in the poten- 
tial correspond to the simple roots Cj — e^+i, where ei, . . . ,eAr denote the 
standard basis elements in M.^ . More generally, the quantum Toda lattice 
associated with a real split semisimple (or reductive) Lie algebra g with 
Cartan subalgebra a has Hamiltonian given by 

where is the Laplacian on a, 11 is a set of simple roots in a* and are 
rational numbers with a particular property [20]. For example, if g = S02Ar+i 
then we can identify a with M^, take 



n = {ei - 62, 62 - 63, ... , CN-l - 6Ar, Cn}, 
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and the corresponding Hamiltonian is given by 



N ^2 

i=l i=l 



The connection between the (generahsed) quantum Toda lattice and the rep- 
resentation theory of the corresponding Lie algebra Q was first observed by 
Kostant [35], who showed that its eigenfunctions can be represented as par- 
ticular matrix elements of infinite-dimensional representations of g. In the 
simplest case when q = 5I2 or QI2, the eigenfunctions are given in terms of 
classical Whittaker functions (actually Macdonald functions). For this rea- 
son, they are often called g- Whittaker functions, or G- Whittaker functions, 
if g = Lie(G). They also arise in the harmonic analysis of automorphic forms 
on Lie groups (see, for example, [8]). There is a spectral decomposition the- 
orem in the general setting due to Semenov-Tian-Shansky [52]. In this paper 
we will only consider the case g = g/jy- However, many of the constructions 
given throughout the paper have Lie-theoretic interpretations and extend to 
the more general setting. This will be indicated where appropriate. 

The eigenfunctions of H have the following integral representation [22, 
28, 19, 26]): 

„ Af-l k 

(4) M^)= / e-^^(^) nn^^'^.- 

Jt(x) 



where T{x) denotes the set of real triangular arrays (T^^j, 1 < i < A; < A^) 
with T/v,i = Xi, 1 < i < N , and 

N / k \ 7V-1 k 

Fx{T) = Afc {Y,Tk,r - Y.^k-i,i ]-Y.Y. (e^'=''-^'=+^-^ + e^'=+i-'+i-^'=.») . 

fe=l \i=l i=l / fe=l i=l 

This integral is defined for A € and has a recursive structure which 
we will now describe. Write H = H^^\ ipx = i^^^^ ■ We will drop these 
superscripts again later, whenever they are unnecessary. For convenience we 
define if(^) = cf/dx"^ and ^p''^\x) = e-^^. Following [19], for TV > 2 and 
9 eC, define a kernel on x M^"^ by 

Qe^\x,y) = exp 

Denote the corresponding integral operator by Qg" ' , defined on a suitable 
class of functions by 
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Then 

V^'' ^Ai,...,Ajv ~ ^Ajv ^Ai,...,Ajv_i' 

and the integral formula (4) can be re-expressed as 

Vx - ^Aat ^Ajv_i ^A2 VXi ■ 

Moreover, as remarked in [19], the following intertwining relation holds: 
This follows from the identity 

which is readily verified. Combining (5) with the intertwining relation (6) 
yields the eigenvalue equation: 

We note the following immediate consequences of the above integral rep- 
resentation. If A G lR^ then ipxix) = i^-xix); if A G lR^ and i/ G R^, then 
|V'A+i/(3;)| < tpuix). It is also known (combining results from [21, 25, 31]) 
that for each x G M^, ipx{x) is an entire function of A G C^. 

The above construction has a representation-theoretic interpretation which 
is described in [19]. It is closely related to the Gauss decomposition and has 
been extended to the other classical Lie algebras in [20]. Encoded in the 
integrand are the defining hyperplanes of the Gelfand-Tsetlin polytope as- 
sociated with the vector x. 

In the present setting (see, for example, [31]), the spectral decomposition 
theorem states that the integral transform 

(7) /(A) = / fix)Mx)dx 

defines an isometry from L2{R'^ , dx) onto L2(iM^, SN{X)dX), where sj\[{X)dX 
is the Sklyanin measure defined by 

(8) -^(^) = 7^^^nr(A,-A.)-^ 
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There is also a Mellin-Barnes type integral formula for -0^ due to Kharchev 
and Lebedev [31, 32, 33] (see also [26]). This is a kind of dual of the above 
integral representation and has a similar recursive structure. For N > 2 and 
z E M, define a kernel on x C^~^ by 

)(A,7) = e^(J^^'-E^')nr(A. - 7,)- 

Then, 

(9) Vr(x) = J Q(f(A,7)<^-^)(x2,...,x^)s^_i(7)d7, 

where the integral is along vertical lines with K7j < KAj for all This 
construction also has a representation-theoretic interpretation which is de- 
scribed in [18]. Gerasimov et al [21] give a clear account of the nature of 
the duality between the two constructions and, in particular, show how this 
duality yields the following identity: for A, z/ E C^, 

(10) / e-''''~'Mx)M^)dx = e^S(^»+^') n r(A, + z/,). 

This is closely related to a Whittaker integral identity which was conjectured 
by Bump [9] and later proved by Stade [55, Theorem 1.1]; there is an exten- 
sive literature on similar and related identities, see for example [10, 27, 54]. 
When N = 2, the eigenfunctions tpx are given by 

M^) = 2exp Q(Ai + A2)(xi + X2)) i^A,-A. (26^^^--^^^) , 

where K^(z) is the Macdonald function. In this case, the Givental's formula 
is equivalent to the integral formula 

Kuiz) = ^^"t^-^exp (-|(t + l/t)) dt, 

the contour integral representation (9) is equivalent to 

K^(z) = — r{s)r{s - u) (-] ds, a > max{5Rz^, 0} 

4'^'' Ja-Loo ^2/ 

and the integral transform defined by (7) is essentially (up to a change of 
variables) the Kontorovich-Lebedev transform. 
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3. The main result. For x,v G M , denote by the probability 
measure on the set F of real triangular arrays {Tk^i)i<i<k<N defined by 



N~l k 



[ fdal = Mx)-' [ /(r)e^''(^) n n 

For i = 1, . . . , N — 1, and continuous rj : (0, oo) — ?• M^, define 

{T^vm = r,{t) + (^log^^e^'+^^^^-^^^^^ds) (e, - e^+i), 

where ei, . . . ,ejv denote the standard basis vectors in M^. Let Hi be the 
identity mapping (IIiTy = rf) and, for 2 < /c < A^— 1, Ilfc = Tio- ■ -oTk-io^k-i- 
Finally, we define 

T = XIat = (7] o • • • o Tn-i) o • • • o (71 o Ti) o Tl. 

The main result of this paper is the following. 

Theorem 3.1 1. If {W{t),t > 0) is a standard Brownian motion in 

with drift v, then {TW{t),t > 0) is a diffusion process in with 
infinitesimal generator given by 

Cu = l^^-' (^^-Jl V'. = + V log • V. 

2. For each t > 0, the conditional law of {{IlkW)i{t), 1 < i < k < N}, 
given {TW{s), s < t; TW{t) = x}, is given by af,. 

3. For each t > 0, the conditional law of W{t), given {TW{s), s < 
t] TW{t) = x}, is given by where 



4- If fJ't denotes the law ofTW{t), then 

Utidx) = exp - - ^ z/ft ipy{x)'dt{x)dx, 

where 

(11) Mx)= I V-A(^)e^»"'*/2^7v(A)dA. 
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It is easy to see that the process {{IlkW)i{t), 1 < i < k < N}, t > 0) 
is Markov. Indeed, setting j = (11^1^) j, it fohows from the construction 
that Z is a Markov process taking values in T which satisfies the system of 
stochastic differential equations: dZi^ = dWi and, for k = 2, . . . , N , 

dZko = dZk-io + (e^^-3-^'--i-2 _ g^fe,2-^fc-i,i) dt 



dZk,k-i = dZk-i,k-i + (e^'^>'--^'=-i''=-i - e.^k,k-i-Zk-i,k-2-^ 
dZk,k = dWk - e^'=-'^-^''-i-'^-i(ii. 

The infinitesimal generator of this process is given by 



where 



l<i<k<N l<i<k<l<N ■ ' l<i<k<N ' 

61,1 (z) = z/i; 

bk,k{^) = ^k- e'^>^-'^-^^>^-\ k = 2,...,N; 

hk,i{z) = e''^'^-''^-^-\ k = 2,...,N; 

hk^i{z) = e^fe.'+i-^fc-i." - (.^k,i-zk-l,^-l ^ l<i<k< N. 

The main content of Theorem 3.1 is the fact that Z^r . is a Markov process, 
with respect to its own filtration. The reason it holds is because 

(12) o = o Au 
where is the Markov operator defined by 

(13) S,/(x) = / f[z)al[dz). 



There is an additional (and non-trivial) technical issue related to the fact 
that these processes start at a particular entrance law coming from '— oo', 
but the intertwining relation (12) lies at the heart of the proof. Actually, the 
proof of Theorem 3.1 given below is based on some intermediate intertwining 
relations which exploit the recursive structure of the quantum Toda lattice 
and the intertwining relation (12) is obtained as a consequence but it should, 
nevertheless, be regarded as the analytic counterpart of Theorem 3.1. As far 
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as we are aware, the intertwining relation (12) and its intermediaries (given 
in section 7 below) have not been previously considered in the literature. 

The operator T was introduced (using a different notation) in the pa- 
per [39], where it was surmised, based on heuristic arguments, that TW 
should be a diffusion process which has the same law as a Brownian mo- 
tion conditioned, in an appropriate sense, on the asymptotic behavior of its 
exponential functionals. In [3] it was observed that such a conditioned Brow- 
nian motion can be defined and moreover, is closely related to the quantum 
Toda lattice, thus providing the impetus for the present work. The above 
notation used to define T follows a more general framework which has been 
developed in the papers [5, 6]. It is shown in [5] that the operators 71 satisfy 
the braid relations, that is, 

71 o 71+1 o 7i = 7i+i o 7i o 7;+!, l<i<N. 
It follows that for each element cr G S jy we can uniquely define 

Ta = Ti^o ■ ■ ■ oTi^ 

where a = + 1) • • • {ip,ip + 1) is any reduced decomposition of a as 

a product of adjacent transpositions. The operator T corresponds to the 
longest element of Sat, that is, 7" = 7^o where 

_ / 1 2 ■■■ N \ 

iv iv-1 ••• I )■ 

The mapping 

Vm ^ imkvUt), l<i<k< N},{Tvis), s < t}) 

is a geometric variant of the RSK (Robinson-Schensted-Knuth) correspon- 
dence. We will explain this connection later and give an interpretation of 
the measure 7g appearing in Theorem 3.1 as a geometric analogue of the 
Duistermaat-Heckman measure associated with the point x. The definition 
of the operator T extends naturally to other Lie algebras, with & j\f replaced 
by the corresponding Weyl group [5, 6]. It is natural to expect the analogue 
of Theorem 3.1 to hold in this more general setting. 

4. The law of the partition function. By Brownian scaling, it is 
easy to see that the processes {Z^{l3),t > 0) and (/3-2(^-i)Z^^(l), t > 0) 

are identical in law, so for convenience we will define Z^^ = Z^(l). The 



DIRECTED POLYMERS AND THE QUANTUM TODA LATTICE 9 

transformation TW is related to the random variable as follows. We 
first note that T satisfies (cf. [5, Lemma 4.6]) 

(14) (_^o)or = ro(-ao), 

where —ao{r]i, . . . , t^at) = (— t/at, . . . , —rji) "Hi denotes the i*'^ coordinate 
of the path rj. It is straightforward to see from the definition of T that 

{TW)N{t) = - log / e-{W,{s^)+W2(suS2)+-+W^(s^-i,t))^^^ ^ ^ ^ 

From the relation (14) we have 

{TW)i{t) = log [ e^N{si)+W^-i{si,S2)+-+W,is^.^,t)^^^ _ _ _ 

-'0<si<-<sjv_i<t 

Thus, if we set W = {Bn, ■■■,Bi) then log Z/^ = (J'W)i{t) and we deduce 
the following. 

Corollary 4.1 The stochastic process (logZf ,t > 0) has the same law as 
the first coordinate of the diffusion process in with infinitesimal generator 

C = ^ij^^Hijo = ^A + VlogVo- V, 
started according to the entrance law 

Ht{dx) = tpQ{x)i9t{x)dx, t > 0, 
where "dt is given by (11). In particular, for n G M, we have 
P(logZf <u)= i^Lt{{x G : xi < n}). 

Note that the relation (14) also implies that the probability measure is 
invariant under the transformation — ctq. Combining Corollary 4.1 with the 
formula (10), we obtain (after shifting the contours in the integral (11) to 
the left in order to apply Fubini's theorem): 

Corollary 4.2 For s > 0, 

Ee-^^^ = Js^ n r(- A0^e5 ^''sN{X)dX, 
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where the integral is along vertical lines with KAj < for all i. 
The probability measure on tM^ with density proportional to 




i>j i<j 



can be interpreted (up to a factor of ltt) as the law, at time 1/t, of the radial 
part of a Brownian motion in the symmetric space of positive definite Her- 
mitian matrices or, equivalently, the law of the eigenvalues of a 'perturbed 
GUE random matrix' A]\f /\/t + Rn /i-, where is, an N x N GUE ran- 
dom matrix and ii at is a diagonal matrix with entries given by the vector 
7r(A'' — 1, A/' — 3, . . . , 1 — A^) (see, for example, [30]). In particular, it is a 
determinantal point process [29] . The above expression for the moment gen- 
erating function of can thus be written as a Fredholm determinant. It 
will be interesting to relate this, in a suitable scaling limit, to the 'crossover 
distributions' recently introduced in the context of KPZ and the stochastic 
heat equation by Sasamoto and Spohn [48, 49, 50, 51] and Amir, Corwin 
and Quastel [1], building on recent work of Tracy and Widom [57, 58, 59, 60] 
on the asymmetric simple exclusion process. See also [12, 13, 14, 15, 45] for 
related recent developments. 

We conclude this section by remarking that the other coordinates of 
TW{t) can also be interpreted as logarithmic partition functions, as fol- 
lows. Define an 'up/right path' in M x Z to be an increasing path which 
either proceeds to the right or jumps up by one unit. For each sequence 
< si < • • • < sn-i < t we can associate an up/right path from 
(0,1) to {t,N) which has jumps between the points {si,i) and {si,i + 1), 
for i = l,...,A^ — 1, and is continuous otherwise. Then we can write 



and the integral is with respect to Lebesgue measure on the Euclidean set of 
all such paths. There is an analogue of Greene's theorem in this context [5] 
(cf. [34]) which yields a similar formula for the other coordinates, namely, 
for each k = 2, . . . , N , 




where 



E{^) = Bi{si) + B2{S2) - B2{si) + ■■■ + BN{t) - Bn{sn-i) 



(TWUt) + ■■■ + {TW)kit) = log / e^('^i)+-+^('^'^)#i . . . 
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where the integral is with respect to Lebesgue measure on the Euchdean set 
of fc-tuples of non-intersecting (disjoint) up/right paths with respective ini- 
tial points (0, 1), . . . , (0, k) and respective end points {t, N—k+1), . . . ,{t, N). 
An interesting property of this formulation is that it extends naturally to 
the continuum setting of KPZ and the stochastic heat equation. 

5. The case N = 2. When N = 2, the eigenfunctions ip^, are given by 

Mx) = 2exp Q(z^i + z^2)(xi K,,_,, (2e(-^--^y^^ . 

In this case, Theorem 3.1 is equivalent to the following theorem of Mat- 
sumoto and Yor [36, 37]. 

Theorem 5.1 1. Let {B[^\t > 0) he a standard one- dimensional Brow- 
nian motion with drift fi, and define 

Zi^^ = l^'^^^'-^i^'ds. 

Then log Z^'') is a diffusion process with infinitesimal generator 

(—logK (e-^)^ — 
2 dx'^ \dx ^ J dx' 

where is the Macdonald function. 

2. The conditional law of b[^\ given {Zs^\ s < t; zj:^^ = z}, is given by 
the generalized inverse Gaussian distribution 

'^K ^{l / z)~^ e^^ exp (- cosh(x)/z) dx. 

3. The law of z[^^ is given by 

P{Z[^'^ G dz) = 2z-^e^/,(t)K^{l/ z)e-^'^'''^dz, 
where 6r{t) is characterized by the Kontorovich-Lebedev transform 

r°° dr 9 , 

2 / Kx{r)er{t)— = e^ '/^, A G iM. 

The above Kontorovich-Lebedev transform can be inverted to obtain 
Orit) = / Kxir)e^ '/^Xsm{7r\)dX. 
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The probability measure Hr^\dt) = lQ{r)~^9r{t)dt is known as the first 
Hartnian- Watson law [24, 37]. It is also characterized by 

e-''^^^^9r{t)dt = I^{r), u > 0, 

where Ix is the modified Bessel function of the first kind. 







6. The zero-temperature limit. By Brownian scaling, we can write 
down a version of Theorem 3.1 for general /? > 0. We will state this in the 
case of zero drift. For continuous rj : (0, oo) — )• M^, define 



{Tfrj){t) = r/(t)+-^log (^/32^*e/^(^»+i(^)-''»(^))(is^ ( 



Ci-ei+i), i = 1, . . . ,iV-l; 



nf = M; nf = 7;^o...o7;lionf_^, 2<k<N; 
Note that 

-^logZf = {T^WUt) - ^log/32. 



Corollary 6.1 1. If W is a standard Brownian motion in M^, then 
1~^W is a diffusion in with generator A/2 + Vlog^o(/3 • ) ' ^• 

2. For each t > 0, the conditional law of {{U'^W)i{t), 1 < i < k < N}, 
given {T^W{s), s < t; T'^W{t) = x}, is given by a^'^if] ■ ). 

3. For each t > 0, the conditional distribution ofW{t), given 
t; T^W{t) = x}, is given by 7^^(/3 • ). 

4- The law ofT^W{t) is given by /i/32i(/3 • ). 

Letting /3 — )■ 00 we recover the multi-dimensional version of Pitman's '2M — 
X' theorem obtained in [7, 43, 40, 5, 6]. For continuous rj : (0, 00) — )• M.^ , 
with 77(0) = 0, define 

{'PiV)i't) = ilit) - mi {r]i{s) - r]i+i{s)){ei - Ci+i), i = 1, . . . , N - 1; 

0<s<t 

Ti = Id.; Ffc =Pio...opfe_iorfe_i, 2<A:<iV; 

•p = Ftv = (-Pi o . . . o Vn-i) o • • • o (Pi o P2) o 7^1- 

By the method of Laplace, as /3 — t- 00, T^W — t- VW uniformly on compact 
intervals and, for each i > and 1 < i < /c < A^, (nfVF)i(t) {rkW)i{t). 
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For 1 < /c < iV, X'' = {(TkW)i, (TkW)k)- By construction, the stochas- 
tic process X(t) = {X^{t), . . . ,X^{t)), t > 0, is Markov and takes values in 
the Gelfand-Tsethn cone 

GTn = {{x\...,x^) eCiX---xCN : x^+i < x'i < x^+\ l<i<k< N-1}, 
where 

Ck = {x £ R'' : xi> ■■■ > Xk}. 

It is a A^(A^ — l)/2-dimensional Brownian motion with singular covariance 
reflected in GT/v via an explicit Skorohod reflection map. But we do not 
need these facts, and refer the reader to the papers [43, 40, 6] for details. 
From the integral formula (4) we have 



„ / N-l k \ 7V-1 k 

tPQ[j3x) = / exp ( - V V (e^fe.'-^fe+i.' +e'^*=+i''+i"^'='') ) TTTTfirfci 
Jr(0x) \ 'H, ~1 / r- , ~_"T ' 



k=l i=l / k=l i=l 

N-l k \ N~l k 



^ ^NiN-i)/2 I exp - ^ ^ + e/5(n+i.+i-n,.)) JJ ,. 

Write = T^j. As /3 — )• oo, if x S Cat, the integrand converges to 1 if 
(x^, . . . , x^) lies in the Gelfand-Tsetlin polytope 

GTNix) = {{x\...,x^) gGTn : x^ = x}, 

and otherwise. It is well-known (for example, by Weyl's dimension formula) 
that the A^(A^ — l)/2-dimensional Euclidean volume of GTi\j{x) is 



n ^- ^(^)' 

\k=l J 



1 

'x] 



where h is given by (2). It follows that, for x € Cat, 
(15) lim /3-^(^-i)/2^o(/3^) = fn ^'1 K 

/3-*-oo \ / 

Similarly, the probability measure crQ^{(] ■ ) converges as /3 — )• oo to the 
uniform probability measure on GT]\j{x). Putting all of this together, letting 
/3 — 7- oo in the statement of Corollary 6.1, we immediately recover parts 1 
and 2 of the following theorem. 

Theorem 6.1 [7, 43, 40, 5, 6] 
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1. If W is a standard Brownian motion in then = VX is a 
Brownian motion conditioned (in the sense of Doob) never to exit Cm- 

2. The conditional law ofli{t), given {X'^{s), s < t; X^{t) = x}, is 
uniform on GTn{x). 

3. The conditional law ofW{t), given {X^{s), s < t; X^{t) = x}, is 
given by the probability measure k^' which is characterized by 



I 



e^^'y^K^'idy) =['[lk\ 



h{x)h{X) 



Part 3 of the above theorem can be deduced from part 2, noting that 
X^ = X]i=i for each \ < k < N . Comparing this with Corol- 
lary 6.1(3) yields the asymptotic formula 

(1.) li. r"'-"'^*v.(^^) - 

This formula can also be seen as consequence of an alternative representation 
of ip\ as an alternating sum of fundamental Whittaker functions [25, 31, 3]. 
The mapping 

Tim ^ iii^kVUt), l<i<k< N},{rv{s), s < t}) 

is a continuous version of the RSK correspondence [40, 5, 6]. The mapping 

Vlo,t] ^ {{{UkvUt), l<i<k< N},{Tv{s), s < t}) 

is a continuous version of the geometric (or 'tropical') RSK introduced by 
Kirillov [34] (see also [38, 4]). The probability measure is the (normalised) 
Duistermaat-Heckman measure associated with the point x. In this setting 
it can be interpreted, via the Harish-Chandra formula, as the conditional 
distribution of the diagonal of a x CUE random matrix given its eigen- 
values X. The probability measure 7o of Theorem 3.1 can thus be inter- 
preted as a geometric analogue of the Duistermaat-Heckman measure. In 
keeping with this analogy, it is natural to record the following analogue of 
the Littlewood- Richardson rule, which follows from Theorem 3.1(3) (c.f. [6, 
Theorem 5.16(ii)]). For s,t>0, define TsW{-) = W{s + ■) - W{s) and 



Gs,t = (T{TW{r), < r < s; {TtsW){u), 0<u<t}. 
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Corollary 6.2 For each x,y £ 



pN 



ipo{x) ipoiy) Jrn ipoiz) 



where 7^'^ is a probability measure on M which can be interpreted, for 
s,t > 0, as the conditional law ofTW{s + t) given Qs,t, TW{s) = x and 
{TT,W){t) = y. 

When N = 2, (17) is equivalent to the formula 

1 /'°° _irf^i'^2_L,„2wti ^, I' zw\ dt 

T' 



Theorem 6.1, in the case = 2, is equivalent to Pitman's celebrated 
'2M — X' theorem [44], which states that, if Xt, t > 0, is a standard one- 
dimensional Brownian motion, then 2maxo<s<f X^ — Xt, t > 0, is a three- 
dimensional Bessel process. Setting W = (i?7v, . . . , i?i) as before, the random 
variable Mi defined by (1) can be written as Mf = (1). Thus, we also 
recover the fact [2, 23] that has the same law as the largest eigenvalue 
of a X GUE random matrix. 

Theorem 6.1 has been generalized to arbitrary finite Coxeter groups in 
the papers [5, 6]. The definition of the operator T also extends naturally 
to other Lie algebras, with ©at replaced by the corresponding Weyl group. 
This is described in [5, 6], where various Lie-theoretic interpretations are 
given. It is natural to expect the analogue of Theorem 3.1 to hold in this 
more general setting. 

7. Intertwining relations. Consider the following extension of the op- 
erator Qg^\ defined on a suitable class of functions / : x M^^^ — )• M 
by 

4'^V(^)= / Qi''\x,y)f{x,y)dy. 
By a straightforward calculation, we obtain 

(18) (i^(^)-Oo7^(^)=7^(^)oc/f), 
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where 



N-l q2 N q2 

i=l "Vi i=l i=l 

d 

OXi 



(9 

_|_ 2 /gj/iv-i-aijv-i _ g2:iV-i-3/jv-2'\ 

aj;Ar_i 

Further integration by parts yields 

(19) (i/(^)-^2)^^W^^(iV)^^W^ 



where 



N-l q2 N q2 

i=l "Vi i=l i=l 

V^yi J dxi 

+ 2 + e"^-?'^ - e"2-^i^ ^ 



9yAr_i y dxN- 

+ 2(9 -e'^'-y^-^) 

^ ^ OXm 

The intertwining relation (19) lies at the heart of this paper. 



8. Proof of Theorem 3.1. We begin by using the intertwining relation 
(19) to prove a Markov functions result. We will then proceed by induction 
to prove a version of Theorem 3.1 for general starting position. The final 
step will be to let the starting position xq — )• — oo (in a sense that will be 
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made precise later). Let v £ M^, and define 



We consider a Markov process {{X{t),Y{t)),i > 0) taking values in x 
defined as follows. The process Y evolves as an autonomous Markov 
process with infinitesimal generator C^^,,^},j^_-i. Let W be standard one- 
dimensional Brownian motion, independent of Y, and define the evolution 
of the process X via the stochastic differential equations 

dXi = dYi + e^^"^^dt 

dX2 = dY2 + (e^3-i2 _ e^2-n) 

dXN^i = dYN^i + (e^^-^^-i - 
dXN = dW+ {uN - e^^-^^-i) dt. 

Then (X, Y) is a Markov process taking values with generator 

Consider the Markov operator A^/^^ defined, for bounded measurable func- 



tions on X ]R(^-i), by 
For X G , define a probability measure on x M^^"^) by 



By (19), we have the intertwining relation 

£f)oAW=AWogf). 
From the theory of Markov functions [47] , we conclude the following. 
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Proposition 8.1 Fix xq,v G M and let (X, y) be a Markov process with 
infinitesimal generator , started with initial law Xf," . Then X is a Markov 
process with infinitesimal generator c[!^^ , started at xq . Moreover, for each 
t >0, the conditional law ofY{t), given {X{s), s < t; X{t) = x}, is given 
by 

The next step is to deduce, by induction, an analogue of Theorem 3.1 for 
general starting position. We construct a Markov process Z taking values in 
r as follows. Let Vl^ be a standard Brownian motion in M.^ with drift ly. The 
evolution of Z is defined recursively by dZi^i = dWi and, for k = 2, . . . , N , 

dZk,i = dZk-1,1 + e^^'^-^^-^'^dt 

dZk,2 = dZk-1^2 + {e^"'-'-^"-^''' - e^^'^-^fe-i.i) dt 

dZk^k-i = dZk-i^k-i + [e^k.k-Zk^i.k-i _ gZfc,fc-i-Zfe-i,fc_2) 
dZk,k = dWk - e^'^'^'^'^-^'^-^dt. 

Proposition 8.2 Fix xq,v € and let Z be the process defined as above 
with initial law a^" . Then Zj\f^. is a Markov process with infinitesimal genera- 
tor C^J^\ started at xq. Moreover, for each t > 0, the conditional law of Z{t), 
given {Zjv,.(s), s < t; Zjv^.(t) = x}, is given by af;, and the intertwining 
relation (12) holds. 

Next we give a formula for the process Z started at Z{Q) = z in terms 
of the driving Brownian motion W . For i = 1, . . . , N — 1, and continuous 
ri : (0,oo) ^ M^, define 

(rSm = r]{t) + log (^e« + 1^ e^"+i(^)-''»(^)ds) (e, - e^+i). 
Fix z gT and, for l<i<k<N — 1, define 

Let n| be the identity map and, for 2 < k < N , 

n| = (r/-^'^o...o7f!r^'=-)on|_i. 

Then, for 1 < i < A; < A^, we can write 

ZkAi) = ^1,1 + (n|VF)i(t). 
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For convenience we will write = 11^ and note that Z^v,. = zi,il + T^W, 
where 1 = (1, 1, ... , 1). Proposition 8.2 can now be restated as follows. 

Proposition 8.3 Fix xq,v G M^. Let W he a standard Brownian motion 
in with drift v and C, a random element of T chosen according to the 
distribution , independent ofW. Then Z^,- = Ci.il + T''W is a Markov 
process with infinitesimal generator c[!^^ , started at xq . Moreover, for each 
t > 0, the conditional law of Z{t), given {Zn,-{s), s < t\ Zj\f^.{t) = x}, is 
given by a^. 

For k = 1, . . . , N, define 



We remark that the vector is half the sum of the positive roots associated 
with the Lie algebra glf.. To complete the proof of Theorem 3.1 we will 
consider the starting position xq = —Mp^, and let M — )■ oo. For this we need 
to understand the asymptotic behavior of ■il)^{—Mp^) and the probability 

AT ^ 

measures ^ as M — )• oo. It was shown by Rietsch [46, Theorem 10.2] 
that the function —Tq{T) on T{x) has a unique critical point T^, which 
is a minimum, and that the Hessian is everywhere totally positive. It is 
straightforward to verify from the critical point equations that 

k N 

k^ '''' - - 

i=l i=l 

Define Si,{T) = J-y{T) — J-q{T) and consider the change of variables 

Tl^, = Tk,i + Mpl l<l<k<N. 
Then we can write 

„ N~l k 

■^r(o) r=i 

It follows, by Laplace's method (see, for example, [17, Theorem 4.14]) that 
the following asymptotic equivalence holds: 

(20) M-Mp"") ~ Ce-^'/^eMe^^^4''\T')) 



k-1 k-1 _ ^ k-l 



as M — )• oo, where C is a constant which is independent of i/. Moreover, 
recalling the above change of variables we see that, in probability, Ck,i — 
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—7- — oo for each l<i<k<N — 1 and Ci,i — ^ 0. It follows by 
the continuous mapping theorem that that Ci,il -'rT'^W converges in law to 
TW, and, for each t > 0, {(n^,VF)j(t), I < i < k < N} converges in law 
to {{nkW)i{t), 1 < i < k < N}. We conclude that TW is a diffusion with 
generator c\!^\ and that the conditional law of {{JlkW)i{t), I < i < k < 
N}, given {TW{s), s < t; TW(t) = x}, is cr^. This proves parts 1 and 2 of 
the theorem. Part 3 of the theorem follows from part 2, noting that for each 
k< N, 

k k-1 
i=l 1=1 

Part 4 follows from part 3 by the spectral decomposition theorem. □ 

Remark 8.1 The asymptotic equivalence (20) is well-known in the case 
N = 2 and can he compared to the full asymptotic expansion obtained in [11] 
in the case N = 3, where it was remarked that the leading term in the 
expansion is independent of the parameter u. 

9. A symmetric version of Proposition 8.2. Proposition 8.2 has a 
'symmetric' analogue which can be regarded as a geometric version of a result 
of Dubedat [16] in the case N = 2, and Warren [61] in the general case. It is 
obtained by applying the intertwining relation (18) rather than (19). In this 
case, we construct a Markov process 5 on P as follows. Let {Wk^i, 1 < i < 
k < N} be a collection of independent standard one-dimensional Brownian 
motions. The evolution of S is defined recursively by dSi,i = dWi^i and, for 
k = 2,...,N, 

dSk,i = dWk,i + e^^-^'^-^^-^dt 

dSk,2 = dWk,2 + (e'^'=-i'2-5fc,2 _ gSfc.2-5fc-i,i) dt 

dSk,k-i = dWk,k-i + [e^k-i,k-i-Sk,k-i _ ^St,k-i-Sk-i,k-2^ 
dSk,k = dWk,k + {i^k - e^^'^-^'^-^-"-^) dt. 

Proposition 9.1 Fix xq, v G and let S be the process defined as above 
with initial law a^" . Then Sn^. is a Markov process with infinitesimal gen- 
erator c[!^\ started at xq. Moreover, for each t > 0, the conditional law of 
S{t), given {S'Ar^.(s), s < t; SN^.(t) = x}, is af,. 

In the case N = 2, with zero drift, we deduce the following. 
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Corollary 9.2 Let Bi, B2 and B^ be independent standard one- dimensional 
Brownian motions. Define 

X{t) = Bi{t)+\og Te^^^^^-^iWds, 

Y{t) = Bsit) - log re^3{^)-B2{.)^g, 
^0 

Then (X + Y)/\/2 is a standard Brownian motion and {X — Y)/\/^ is a 
diffusion process (independent of {X -\-Y)/^/2) with infinitesimal generator 

2 dx'^ \dx J dx 

Acknowledgements. The author would like to thank the referees for careful 
reading and helpful suggestions which have led to an improved version of the 
paper. Thanks also to Philippe Biane, Philippe Bougerol, Ivan Corwin, Timo 
Seppalainen and Herbert Spohn for helpful and stimulating discussions, and 
Dan Bump for pointers to the literature on Whittaker integral identities. 

REFERENCES 

[1] G. Amir, I. Corwin and J. Quastel. Probability distribution of the free energy of the 
continuum directed random polymer in 1+1 dimensions. Comm. Pure Appl. Math., 
to appear. 

[2] Y. Baryshnikov. GUEs and queues. Probab. Th. Rel. Fields 119 (2001) 256-274. 

[3] F. Baudoin and N. O'Connell. Exponential functionals of Brownian motion and class 

one Whittaker functions. Ann. Inst. H. Poincare Probab. Statist., to appear. 
[4] A. Berenstein and D. Kazhdan. Geometric crystals and unipotent crystals. GAFA 

2000, Tel Aviv, 1999, Geom. Funct. Anal. Special Volume, Part I (2000). 
[5] Ph. Biane, Ph. Bougerol and N. O'Connell. Littelmann paths and Brownian paths. 

Duke Math. J. 130 (2005) 127-167. 
[6] Ph. Biane, Ph. Bougerol and N. O'Connell. Continuous crystals and Duistermaat- 

Heckman measure for Coxeter groups. Adv. Math. 221 (2009) 1522-1583. 
[7] Ph. Bougerol and Th. Jeulin. Paths in Weyl chambers and random matrices. Probab. 

Th. Rel. Fields 124 (2002) 517-543. 
[8] D. Bump. Automorphic forms on GL(3,M.). Lecture Notes in Mathematics, 1083. 

Springer- Verlag, Berlin, 1984. 
[9] D. Bump. The Rankin Selberg method: a survey, in Number Theory, Trace Formulas, 

and Discrete Groups (K. E. Aubert, E. Bombieri and D. Goldfeld, eds.). Academic 

Press, New York, 1989. 
[10] D. Bump and S. Friedberg. The exterior square automorphic L-functions on GL{n). 

Festschrift in Honor of Piatetski-Shapiro, Part II, Weizmann, Jerusalem, 1990, pp. 

4765. 



22 N. O'CONNELL. 

[11] D. Bump and J. Huntley. Unramified Whittaker functions for GI/(3,R). J. Anal. 

Math. 65 (1995) 19-44. 
[12] P. Calabrese, P. Le Doussal and A. Rosso. Free-energy distribution of the directed 

polymer at high temperature. EPL 90 (2010) 20002. 
[13] V. Dotsenko. Replica Bethe ansatz derivation of the Tracy- Widom distribution of the 

free energy fluctuations in one-dimensional directed polymers. J. Stat. Mech. (2010) 

P07010. 

[14] V. Dotsenko. Bethe ansatz derivation of the Tracy- Widom distribution for one- di- 
mensional directed polymers. EPL 90 (2010) 20003. 

[15] V. Dotsenko and B. Klumov. Bethe ansatz solution for one-dimensional directed 
polymers in random media. J. Stat. Mech. (2010) P03022. 

[16] J. Dubedat. Reflected Brownian motions, intertwining relations and crossing proba- 
bilities. Ann. Inst. H. Pomcare Probab. Statist. 40 (2004) 539-552. 

[17] M. Evans and T. Swartz. Approximating Integrals via Monte Carlo and Deterministic 
Methods. Oxford University Press, 2000. 

[18] A. Gerasimov, S. Kharchev and D. Lebedev. Representation Theory and the Quantum 
Inverse Scattering Method: The Open Toda Chain and the Hyperbolic Sutherland 
Model. Int. Math. Res. Notices 2004, 823-854. 

[19] A. Gerasimov, S. Kharchev, D. Lebedev and S. Oblezin. On a Gauss-Givental rep- 
resentation of quantum Toda chain wave equation. Int. Math. Res. Notices 2006, 
1-23. 

[20] A. Gerasimov, D. Lebedev and S. Oblezin. New Integral Representations of Whittaker 

Functions for Classical Lie Groups. arXiv:0705.2886. 
[21] A. Gerasimov, D. Lebedev and S. Oblezin. Baxter Operator and Archimedean Hecke 

Algebra. Commun. Math. Phys. 284 (2008) 867-896. 
[22] A. Givental. Stationary phase integrals, quantum Toda lattices, flag manifolds and 

the mirror conjecture. Topics in Singularity Theory, AMS Transl. Ser. 2, vol. 180, 

AMS, Rhode Island (1997) 103-115. 
[23] J. Gravner, C.A. Tracy and H. Widom. Limit theorems for height fluctuations in a 

class of discrete space and time growth models. J. Stat. Phys. 102 (2001) 1085-1132. 
[24] P. Hartman and G. Watson. "Normal" distribution functions on spheres and the 

modified Bessel functions. Ann. Probab. 2 (1974) 593-607. 
[25] M. Hashizume. Whittaker functions on semisimple Lie groups. Hiroshima Math. J. 

12 (1982), no. 2, 259-293. 
[26] T. Ishii and E. Stade. New formulas for Whittaker functions on GL{n,'M.). J. Fun. 

Anal 244 (2007) 289-314. 
[27] H. Jacquet. Integral representations of Whittaker functions. In: Contributions to Au- 

tomorphic Forms, Geometry and Number Theory (Shalika volume), eds H. Hida, D. 

Ramakrishnan and F. Shahidi. Johns Hopkins University Press, Baltimore, MD, 2004. 
[28] D. Joe and B. Kim. Equivariant mirrors and the Virasoro conjecture for flag mani- 
folds. Int. Math. Res. Notices 2003, no. 15, 859-882. 
[29] K. Johansson. Determinantal processes with number variance saturation. Comm. 

Math. Phys. 252 (2004) 111-148. 
[30] L. Jones and N. O'Connell. Weyl Chambers, Symmetric Spaces and Number Variance 

Saturation. ALEA Lat. Am. J. Probab. Math. Stat. 2 (2006) 91-118. 
[31] S. Kharchev and D. Lebedev. Integral representations for the eigenfunctions of a 

quantum periodic Toda chain. Lett. Math. Phys. 50 (1999), 53-77. 



DIRECTED POLYMERS AND THE QUANTUM TODA LATTICE 



23 



[32] S. Kharchev and D. Lebedev. Eigenfunctions of GL{N,'R) Toda chain: the Mellin- 
Barnes representation. JETP Lett. 71 (2000) 235-238. 

[33] S. Kharchev and D. Lebedev. Integral representations for the eigenfunctions of quan- 
tum open and periodic Toda chains from the QISM formalism. J. Phys. A 34 (2001) 
2247-2258. 

[34] A. N. Kirillov. Introduction to tropical combinatorics. Physics and Combinatorics. 

Proc. Nagoya 2000 2nd Intemat. Workshop (A. N. Kirillov and N. Liskova, eds.), 

World Scientific, Singapore, 2001, pp. 82-150. 
[35] B. Kostant. Quantisation and representation theory. In: Representation Theory of Lie 

Groups, Proc. SRC/LMS Research Symposium, Oxford 1977, LMS Lecture Notes 34, 

Cambridge University Press, 1977, pp. 287-316. 
[36] H. Matsumoto and M. Yor. A version of Pitman's 2M — X theorem for geometric 

Brownian motions. C. R. Acad. Sci. Pans 328 (1999) 1067-1074. 
[37] H. Matsumoto and M. Yor. An analogue of Pitman's 2M — X theorem for exponential 

Wiener functionals, part I: a time inversion approach. Nagoya Math. J. 159 (2000) 

125-166. 

[38] M. Noumi and Y. Yamada. Tropical Robinson-Schensted-Knuth correspondence and 

birational Weyl group actions. Adv. Stud. Pure Math. 40 (2004) 371442. 
[39] N. O'Connell. Random matrices, non-colliding processes and queues. Semmaire de 

Probabilites, XXXVI, 165-182, Lecture Notes in Math., 1801, Springer, 2003. 
[40] N. O'Connell. A path-transformation for random walks and the Robinson- Schensted 

correspondence. Trans. Amer. Math. Soc. 355 (2003) 3669-3697. 
[41] N. O'Connell and J. Moriarty. On the free energy of a directed polymer in a Brownian 

environment. Markov Process. Related Fields 13 (2007) 251-266. 
[42] N. O'Connell and M. Yor. Brownian analogues of Burke's theorem. Stock. Process. 

Appl. 96 (2001) 285-304. 
[43] N. O'Connell and M. Yor. A representation for non-colliding random walks. Elect. 

Comm. Probab. 7 (2002). 
[44] J.W. Pitman. One-dimensional Brownian motion and the three-dimensional Bessel 

process. Adv. Appl. Probab. 7 (1975) 511-526. 
[45] S. Prolhac and H. Spohn. Two-point generating function of the free energy for a 

directed polymer in a random medium. arXiv:1011.4014. 
[46] K. Rietsch. A mirror construction for the totally nonnegative part of the Peterson 

variety. Nagoya Math. J. 183 (2006) 105-142. 
[47] L. C. G. Rogers and J. Pitman. Markov functions. Ann. Probab. 9 (1981) 573-582. 
[48] T. Sasamoto and H. Spohn. Exact height distributions for the KPZ equation with 

narrow wedge initial condition. Nucl. Phys. B 834 (2010) 523-542. 
[49] T. Sasamoto and H. Spohn. One-dimensional Kardar-Parisi-Zhang equation: An exact 

solution and its universality. Phys. Rev. Lett. 104, 230602 (2010). 
[50] T. Sasamoto and H. Spohn. The crossover regime for the weakly asymmetric simple 

exclusion process. J. Stat. Phys. 140 (2010) 209-231. 
[51] T. Sasamoto and H. Spohn. The l-fl-dimensional Kardar-Parisi-Zhang equation and 

its universality class. Proceedings StatPhys 24; Journal of Statistical Mechanics, on- 
line (2010). 

[52] M. Semenov-Tian-Shansky. Quantisation of open Toda lattices. In: Dynamical sys- 
tems VIL Integrable systems, nonholonomic dynamical systems. Edited by V. I. 



24 



N. O'CONNELL. 



Arnol'd and S. P. Novikov. Encyclopaedia of Mathematical Sciences, 16. Springer- 
Verlag, 1994. 

[53] T. Seppalainen and B. Valko. Bounds for scaling exponents for a 1+1 dimensional 
directed polymer in a Brownian environment. arXiv: 1006.4864 

[54] E. Stade. Mellin transforms of GL{n, K) Whittaker functions. Amer. J. Math. 123 
(2001) 121-161. 

[55] E. Stade. Archimedean L-factors on GL{n) x GL{n) and generalized Barnes integrals. 

Israel J. Math. 127 (2002) 201-219. 
[56] C.A. Tracy and H. Widom. Level-spacing distributions and the Airy kernel. Comm. 

Math. Phys. 159 (1994) 151-174. 
[57] C.A. Tracy and H. Widom. Integral formulas for the asymmetric simple exclusion 

process. Comm. Math. Phys. 279 (2008) 815-844. 
[58] C.A. Tracy and H. Widom. A Fredholm determinant representation in ASEP. J. Stat. 

Phys. 132 (2008) 291-300. 
[59] C.A. Tracy and H. Widom. Asymptotics in ASEP with step initial condition. Comm. 

Math. Phys. 290 (2009) 129-154. 
[60] C.A. Tracy and H. Widom. Formulas for joint probabilities for the asymmetric simple 

exclusion process. J. Math. Phys. 51 (2010) 063302. 
[61] J. Warren. Dyson's Brownian motions, intertwining and interlacing. Electron. J. 

Prohah. 12 (2007) 573-590. 



Mathematics Institute, University of Warwick, Coventry CV4 7AL, UK, 
E-MAIL: n.m.o-connellOwarwick. ac.uk 



